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We investigate the survival rate of an initial momentum anisotropy (u 2 "’)i not spatial anisotropy, 
to the final state in a multi-phase transport (AMPT) model in Au-|-Au collisions at ,/snn=2Q0 GeV. 
It is found that both the final-state parton and charged hadron V 2 show a linear dependence versus 
respect to the participant plane (PP). It is found that the slope of this linear depen¬ 
dence (referred to as the survive rate) increases with transverse momentum (pr), reaching ~100% 
at pt^2.5 GeV/c for both parton and charged hadron. The survival rate decreases with collision 
centrality and energy, indicating decreasing survival rate with increasing interactions. It is further 
found that a u^^'IRnd} with respect to a random direction does not survive in W 2 {PP} but in the 
two-particle cumulant V2{2}. The dependence of V2{2} on U 2 ’^*{Rnd} is quadratic rather than linear. 


I. Introduction 

A new state of matter, the strongly coupled quark 
gluon plasma (sQGP), is created in relativistic heavy ion 
collisions [1-4]. One of the most important evidence is 
the measured large elliptic flow in non-central heavy ion 
collisions, believed to stem out of final state interactions 
in the anisotropic overlap zone [5]. The measured ellip¬ 
tic flow is so large that it is compatiable with hydrody¬ 
namic calculations with minimal shear viscosity to en¬ 
tropy density ratio (rj/s), indicating maximal final-state 
interactions [6, 7]. 

Present hydrodynamic calculations start from an ini¬ 
tial condition of isotropic momentum distribution. It 
has been argued, however, that the initial momentum 
anisotropy may not be zero in relativistic heavy ion colli¬ 
sions. For example, it is suggested that the wave function 
is asymmetric in momentum space due to Heisenberg un¬ 
certainty principle because of the spatial anisotropic over¬ 
lap [8] . In classical Yang-Mills dynamics it is found that 
initial momentum anisotropy can arise from the event-by¬ 
event breaking of rotational invariance in local domains 
whose size is determined by the saturation scale [9]. Ini¬ 
tial flow in classical Yang-Mills field can also develop from 
the non-abelian generalization of Gauss’ Law and Am¬ 
pere’s and Faraday’s Laws [10]. In proton-proton colli¬ 
sions color reconnection can produce initial flow-like cor¬ 
relations [11] and it may be relevant for heavy ion col¬ 
lisions as well. If there indeed exist initial momentum 
anisotropies and these initial anisotropies can partially 
survive to the final state, then the comparison of data 
to hydrodynamics without initial momentum anisotropy 
would not be reliable to extract transport properties of 
the sQGP, such as the 77 /s. In this paper, we investigate 
the survival rate of an input initial momentum anisotropy 
using a parton transport model. 
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II. Analysis 

We employ A Multi-Phase Transport (AMPT) model 
with string melting and 3 mb parton cross section [12, 13]. 
This model can describe well the measured particle ra¬ 
pidity distributions, transverse momentum spectra, and 
elliptic flow [14]. AMPT consists of four main parts: the 
initial condition, parton-parton interactions, hadroniza- 
tion, and hadronic scatterings. The initial condition is 
taken from the HIJING model [15]. It uses Glauber 
nuclear geometry to model the spatial and momentum 
information of minijet partons from hard processes and 
strings from soft processes. The interactions of partons 
are treated by the ZPG parton cascade model [16]. Af¬ 
ter parton interactions cease, a combined coalescence and 
string fragmentation model is used for the hadronization 
of partons. Finally, the ART model is used to describe 
the elastic and inelastic scatterings of hadrons [17]. 

Elliptic flow can be quantified by V 2 , the second har¬ 
monic Fourier coefficient of the particle azimuthal distri¬ 
bution in momentum space [18], 

dN/d(j) oc 1 -b 2v2{PP} cos 2((/ - 4'2{PP}). (1) 

In AMPT, the initial parton azimuthal distribution is 
isotropic: 

dN 

—— = constant. (2) 

We can artificially create a momentum anisotropy by 
“squeezing” particles towards a particular plane. We first 
choose this plane to be the participant plane (PP) of the 
initial partons in configuration space. The azimuthal an¬ 
gle of the participant plane is given by 

atan2((r2„^ sin2(/i„i), cos 2 </™)) + ,o^ 

where rini and (fini are polar coordinate position. Mathe¬ 
matically we change each parton’s initial azimuthal angle 
(fini into by: 


4^1711 — /’irii + 


( 4 ) 
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such that 


dN 


= l + 2ur{PP}cos2(C.-'I^r)- (5) 


In order to achieve an initial anisotropy u^HPP} with 
respect to d' 2 ^, one applies 


5 = - 


ur{PP}sin2(0,„, 

1 + 2ur{PP} cos2(<^,™ - 


( 6 ) 


Second, we choose this particular plane to be a random 
azimuthal direction participant plane along 

which hydrodynamic collective flow develops. We denote 
this initial anisotropy as U 2 ”*{Rnd}. Same as Eq. 6 , on 
applies 


u^"HRnd}sin2((/>„, 

1 + 2u^™{Rnd} cos 2((/>,„, - ' 


In these operations, only the parton’s azimuthal angle 
is altered, no other changes. The event now has an ini¬ 
tial anisotropy (u^'lPP} or ?; 2 "*{Rnd}). The event then 
evolves as modeled by AMPT. In this analysis we have 
used a given u™* in each event, independent of the parton 
Pt- 

We analyze the momentum anisotropies of the final- 
state partons (i.e. after parton interactions cease and be¬ 
fore hadronization) and the final-state hadrons by the 
Fourier coefficients [19]: 

U2{PP} = (cos2(<^-vI/PP)), (8) 


where </> is the particle (parton or hadron) azimuthal an¬ 
gle. 

Experimentally, however, the participant plane is in¬ 
accessible. The momentum anisotropy is often analyzed 
by final-state two-particle correlations. In absence of in¬ 
trinsic particle correlations (nonflow), the final-state two- 
particle correlations are caused by each particle’s corre¬ 
lation to the common participant plane (i.e. flow corre¬ 
lations) of Eq. 1. In such a case, the two-particle density 
is given by 

d^N/dA4> =1-1- 2u2{2}^ cos A^, (9) 


where Acj) = ipi — (j >2 is the azimuthal angle difference 
between the two particles. The momentum anisotropy 
can thereby be calculated by 

U2{2} = v'((cos2(0i -(/>2))), (10) 

without requiring the participant plane [20]. We use the 
Q-cumulant method [21] to calculate (cos2((/)i — (j) 2 )) in 
each event, by 

Here the second-harmonic flow vector Q 2 is given by 

M 

Q2=^e*2^% (12) 

i 


and M is the number of particles used in the Q-cumulant. 
In the averaging of (cos 2(0i — 02 ) over events in Eq. 10, 
we have applied the weight of M{M — 1). For differen¬ 
tial flow as a function of pt, we compute the correlation 
between a particle of interest at pt (azimuth 0) and a 
reference particle (azimuth 0ref), ((cos 2(0 — 0i.ef))). The 
anisotropy of the particles of interest is then given by 

^ (cos 2 ( 0 - 0ref)) 

=-<{ 2 )- 

where the reference particle V2‘^^{2} is given by Eq. 10, 
where 0i and 02 are both of reference particles. We again 
use the Q-cumulant method to compute (cos 2(0 —0i.ef)), 
by 

= (14) 

where p 2 is the second-harmonic flow vector (Eq. 12) of 
the particles of interest, nip is the number of the particles 
of interest, and is the number of particles in the over¬ 
lap of the two groups (particles of interest and reference 
particles). 


III. Results 

Figure 1 shows the final-state charged hadron U 2 {PP} 
at mid-rapidity (jTyj < 1 ) versus transverse momen¬ 
tum (pt) for fixed impact parameter (&=8 fm) Au-|-Au 
collisions at y/SNN=‘2'00 GeV. The black circles show 
the default result without initial u^'jPP}. The red 
squares show the result with U 2 "®{PP}= 8 %. With ini¬ 
tial U 2 ™{PP}, the final U 2 {PP} is larger, and the effect is 
more significant at higher py. This is consistent with the 
finding in Ref. [22] that higher pT partons suffer fewer 
collisions on average and thus retain larger fraction of 
their initial anisotropy. 

Figure 2 shows the parton and hadron U 2 {PP} versus 
initial parton ^“‘{PP}, for 1.5< pr <2 GeV/c as an ex¬ 
ample. We fit the results with the two-parameter linear 
function, U 2 {PP} = r x U 2 "*{PP} -l-U 2 ( 0 ). The fitting pa¬ 
rameter ^ 2 ( 0 ) corresponds to the result without an initial 
^mzjpp} ^j^g gigpg parameter r to quantify the 

survival rate of an input initial u™®{PP}. We show in 
Fig. 3 the survival rate as a function of pr for &=8 fm 
Au-|-Au collisions. The survival rate increases with py. 
This is because partons with lower px suffer on average 
more collisions before freezeout, which tend to wash out 
the initial V 2 - Meanwhile, at higher pr, the survival rate 
is larger as they suffer fewer collisions. With zero colli¬ 
sion, the initial U 2 "*{PP} will be intact and the survival 
rate would be 100%. It is interesting to note that the 
survival rate at high px can be larger than unity, at least 
for the partons in the highest px bin studied in this anal¬ 
ysis. This probably suggests some non-linear effect in V 2 
at high Px that the initial V 2 enhances the final-state de¬ 
velopment of collective anisotropic flow. It is worthwhile 
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FIG. 1. (color online) Mid-rapidity (|?;| < 1) hadrons n 2 {PP} 
versus pt with U 2 "*{PP} = 0% and U 2 "*{PP} = 8% for b—8 fm 
Au+Au collisions at -^snn=‘200 GeV by AMPT (string melt¬ 
ing)- 


FIG. 3. (color online) The survival rate of initial momentum 
anisotropy to final-state partons and charged hadrons at mid¬ 
rapidity {\rj < 1|) as a function of pt for b—8 fm Au+Au 
collisions at ^snn=‘2Q0 GeV by AMPT (string melting). 


to note that the survival rate of the initial ^“^{PP} to 
the final-state hadron, at a given px value, is smaller than 
that to the final-state partons before hadronization. This 
is presumably due to the facts that partons cluster into 
hadrons at higher px and that hadrons rescatter after 
hadronization. 
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FIG. 2. (color online) Mid-rapidity (|p| < 1) final-state par- 
tons and hadrons U 2 {PP} at 1.5 < pt < 2 GeV/c versus 
^_g Au+Au collisions at ^sjvjv=200 GeV 
by AMPT (string melting). 


We show the beam energy dependence of the sur¬ 
vival rate in Fig. 4 (a) and the centrality dependence in 
Fig. 4 (b). The parton and haron px are integrated over 
the range of 0< px <2 GeV/c. Note V 2 is generally not 
zero in 6=0 fm collisions because of event-by-event fluctu¬ 
ations in the initial overlap geometry. The survival rate 
decreases with increasing collision energy and increasing 
centrality for both partons and hadrons. Higher energy 
collisions and/or more central collisions correspond to 
stronger interactions which reduce the survival rate of 


the initial anisotropy. 

The interpretation of our results is relatively straight¬ 
forward. The survival rate depends on the final-state 
activity. The more the partons (and hadrons) interact, 
the smaller the survival rate of the initial momentum 
anisotropy. Because in AMPT the collision system has 
relatively low opacity [22], the survival rate is apprecia¬ 
ble. With large opacity, the survival rate should be min¬ 
imal. It will be interesting to repeat our study using 
hydrodynamics, starting from a non-zero initial momen¬ 
tum anisotropy. 
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FIG. 4. (color online) Survival rate as a function of (a) 
collision energy and (b) impact parameter for mid-rapidity 
(I 77 I < 1 ) final-state partons and charged hadrons in 0< px <2 
GeV/c for Au+Au collisions by AMPT (string melting). 


We have so far used the participant plane from the 
model to calculate hnal-state particle W 2 {PP} to demon¬ 
strate that initial momentum anisotropy can survive 
to the final state. As mentioned earlier, participant 
plane cannot be accessed in experiment, and instead two- 
particle correlations are used to measure the momentum 
anisotropy, U2{2}. We have verified that U2{2}, while 
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slightly larger than r' 2 {PP} because of nonhow effects in 
AMPT, gives the same qualitative conclusion. 

In the above we have generated v^'^{PP} with respect 
to the participant plane. Several physics mechanisms [9- 
11] suggest that the initial momentum anisotropy is inde¬ 
pendent of the conhguration space anisotropy of the col¬ 
lision, but rather random. To study the survival rate of 
such initial momentum anisotropies, we repeat our anal¬ 
ysis but with n 2 "*{Rnd} with respect to a random plane 
by Eq. 7. 

Figure 5 (a) shows the U 2 {PP} as a function of 
Pt with various u^^jRnd}. As seen from the hgure, 
?;™®{Rnd} does not have major effect on z; 2 {PP}. In 
other words, ri^^iRnd} does not survive in n 2 {PP}. This 
is understandable because the initial correlations due to 
^““{Rnd} are averaged out to zero in W 2 {PP}- Fig¬ 
ure 5 (b) shows the V2{‘2} as a function oipr with various 
n 2 "*{Rnd}. It is found that ti 2 ™{Rnd} has noticeable ef¬ 
fect on the hnal V 2 { 2 }. 
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FIG. 5. (color online) Mid-rapidity (|?;| < 1) hadrons (a) 
r' 2 {PP} and (b) V 2 { 2 } versus pr with U 2 "®{Rnd}= 0 %, 8 % 
and 20% for 6=8 fm Au-|-Au collisions at ^sjviv=200 GeV by 
AMPT (string melting). 


sus i'^ Pt ranges. The V2{‘2} is slightly 

larger than z; 2 {PP} with no r;™*, and this is because 
of nonhow effects in AMPT which are present in W2{2} 
but not in 'y 2 {PP}- The horizontal dashed line (going 
through the hollow point at r'™®{Rnd}=0) would indicate 
vanishing effect of ^^'{Rnd}. As can be seen from the 
hgure, the f 2 {PP} data points approximately he along 
the dashed line. Some minor, but noticeable deviations 
from the horizontal line may be due to the fact that the 
AMPT evolution will not be the same once an initial 
?;™®{Rnd} is introduced in an event. A large z; 2 ”*{Rnd} 
may drive the collision system to expand more efficiently 
in a direction not coinciding with the participant plane 
but slightly deviating from it. This could cause a smaller 
r; 2 {PP}, which seems to be the case for large f^jRnd}. 

The n2{2}, on the other hand, increases with 



FIG. 6 . (color online) Mid-rapidity (I 77 I < 1) hadrons V2{2} 
and U 2 {PP} versus U 2 "*{Rnd} for b =8 fm Au-|-Au collisions 
at y/SNN=200 GeV by AMPT (string melting). The solid 
lines are htting results. The dashed lines are horizontal and 
through the point of ii 2 "'{Rnd}= 0 . 


In order to study the survival rate of ^^"'“{Rnd}, we 
show in Fig. 6 the mid-rapidity (Iryj < 1) charged hadron 
'C 2 { 2 } (solid markers) and t' 2 {PP} (hollow markers) ver- 


served with n 2 "*{PP}. This can be understood as follows. 
Suppose the Hnal single particle azimuthal distribution is 
described by 


c?iV/d(/) oc 1 -I-2i;2{PP}(0) cos(0 —'I '2 )-I- 2r x'u^^lRnd} cos2((/) — ), 


(15) 


where ?;2{PP}(0) is the “hydrodynamic” flow from normal AMPT, and rx?;2"*{Rnd} is the surviving initial momentum 
anisotropy along 'I'2{Rnd} with a rate r. The two-particle correlations would be 

d'^N/dAcj) oc J J -b 2r;2{PP}(0) cos((?ii - T^^) -b 2r x r;™{Rnd} cos2((/.i - 4'^"'^)] 

X [1 -b 2r;2{PP}(0) cos(<^2 - ^'2^) -b 2r x n^“{Rnd} cos2(((.2 - 

oc 1-b 2ti2{PP}(0)^ cos2A^-b x 2n“*{Rnd}^ cos 2A0. (16) 

I- 

Thus, comparing to Eq. 9 , we have The relationship is not linear, but quadratic. If the initial 

V2{2} = y^U 2{PP}(0)2 -b r2 X ?;i™{Rnd}2. (17) 
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ti™* was with respect to the participant plane, as we have 
studied previously, then the final V 2 would be 

V2{2} ^ U 2 {PP} = U2{PP}(0) + r X wr{PP}, (18) 

as we have used in the linear fit in Fig. 2. We thus fit the 
V2{2} vs. ?; 2 "*{Rnd} data in Fig. 6 by a quadratic func¬ 
tion of the form in Eq. 17. The fits are shown by the solid 
curves in Fig. 6. Our results show that ^““{Rnd} will 
survive to final state. This is because the ?; 2 "*{Rnd} in¬ 
troduces an initial particle correlation, which survives to 
the end when the particles interact minimally in the final 
state. Therefore, the initial momentum anisotropy, even 
with respect to a random plane, will contribute to the 
final anisotropy measurement by particle correlations. 


IV. Summary 

We have studied to what extent an input initial mo¬ 
mentum anisotropy survives to the final state using 
the AMPT model. It is found that the final momen¬ 
tum anisotropy shows a linear dependence on the ini¬ 
tial momentum anisotropy relative to the configuration- 
space participant plane (PP). The slope of this lin¬ 


ear dependence (survival rate) increases approximately 
linearly with px, with ^100% survival rate at high 
Pt^2.5 GeV/c. The survival rate decreases with increas¬ 
ing centrality and increasing beam energy. An initial mo¬ 
mentum anisotropy relative to a random plane does not 
survive to the final-state ?; 2 {PP} with respect to the par¬ 
ticipant plane, but survives to the two-particle cumulant 
V2{2}. The V2{2} survival rate is quadratic in n^'jRnd} 
with appreciable magnitude when ri 2 ”*{Rnd} is sizeable. 
Our results can be understood as that the survival rate 
decreases with increasing final-state interactions. 
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